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CONDITIONING SUPER-BROWNIAN MOTION ON ITS BOUNDARY 
STATISTICS AND FRAGMENTATION 
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We condition super-Brownian motion on "boundary statistics" of 
the exit measure Xd from a bounded domain D. These are random 



o 



variables defined on an auxiliary probability space generated by sam- 
pling from the exit measure Xd- Two particular examples are: condi- 
tioning on a Poisson random measure with intensity fiXo', and con- 
ditioning on Xd itself. We find the conditional laws as /i-transforms 
of the original SBM law using Dynkin's formulation of X-harmonic 
functions. We give explicit expression for the (extended) X-harmonic 
functions considered. We also obtain explicit constructions of these 
conditional laws in terms of branching particle systems. For exam- 
ple, we give a fragmentation system description of the law of SBM 
conditioned on Xd = v, in terms of a particle system, called the 
. backbone. Each particle in the backbone is labeled by a measure v, 

0^ ' representing its descendants' total contribution to the exit measure. 

The particle's spatial motion is an /i-transform of Brownian motion, 
where h depends on v. At the particle's death two new particles are 



' born and ii is passed to the newborns by fragmentation. 

1. Introduction. Studying conditioned Markov processes is a kind of inverse problem - given 
information about how the process ends up, one tries to infer how it got there, at least in terms of 
t**- ■ probabilities. In the context of Brownian motion and finite dimensional Markov processes one can 

make very explicit calculations, starting with the work of J. L. Doob (1959). Attempts to make 
\ similar calculations for super Brownian motion are more recent. These studies typically aim to 

■ recover the conditional law of a superprocess as the law of a distinct probabilistic object. Several 

authors have succeeded in coming up with such descriptions for certain conditionings, and produced 
models with remarkably rich structure. The first of these models was the immortal particle system 
of Evans and Perkins (1990) and Evans (1993), for super-Brownian motion in M d conditioned on 
survival. In this model, an immortal particle moves according to a Brownian motion, and throws 
off mass at a uniform rate, and then this mass evolves in the space as an unconditioned super- 
Brownian motion. Following Evans and Perkins, Salisbury and Verzani (1999) considered a super 
Brownian motion X in a domain D, with conditioning based on the exit measure Xo from D. 
More specifically, they conditioned X on the event that the support of Xo contains certain points 
z±,. . . ,Zk, and recovered the resulting conditional law in terms of a branching backbone system. 
The branching backbone is a random tree with k leaves reaching the points z\ . . . z^. Similar to 
Evans and Perkins' model, there is mass uniformly created along the branching backbone which 
follows the law of an unconditioned super-process independent of the points Z\, . . . , z^. Giving such 
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an explicit characterization of a conditioned process is an interesting problem from a probabilistic 
modeling point of view. For example, in population dynamics, one can view it as an analogue of 
a host of biological problems in which one has information about the state of the population at 
certain times or locations, and one wishes to infer the genealogical structure of the populations 
of the ancestors (e.g. the "out of Africa" problem of human origins). For explicit representations 
of other related conditioned processes see Roelly-Coppoletta and Rouault (1989), Overbeck (1993, 
1994), Etheridge (1993), or Etheridge and Williams (2003). 

It turns out that there is more to the conditioning problem than described above. A conditioned 
process represents a special case of a Girsanov transformation, or a martingale change of measure. 
For concreteness, let us consider the following example: Let £ be Brownian motion in a domain 
E. Let te be the exit time from E. We compute the conditional law, U z of £ given £ TB = z, by a 
martingale change of measure from U x , the law of £. This martingale change of measure is given in 
terms of a certain harmonic function h z (-). More precisely, for any domain D compactly contained in 
E, and any Y measurable with respect to J~ TD , we have H%,(Y) = H x (Yh z (^ TD )/h z (x)). A Girsanov 
transformation defined in terms of a harmonic function is called an /i-transform, and typically 
conditional laws of Markov processes are formulated as /i-transforms of their original laws. This 
classical relationship between harmonic functions and conditioning a Markov process leads to an 
elegant probabilistic formulation of the Martin boundary theory for elliptic differential operators. 

In the context of super-processes, the analogue of harmonic functions are X-harmonic functions. 
Following the definition of Dynkin (2002), let us consider a super-Brownian motion X = (Xjj, P^), 
a family of random measures (exit measures) and their associated probability laws where D is an 
open subset of a given domain E in M. d and fi is a finite measure on E. Write D <s E if D is 
open and its closure is a compact subset of E. A non-negative function H is X -harmonic if for any 
D E and any finite measure /i with support in D, 

(1) P,(H(X D )) = H(v). 

Note that this property resembles the mean value property of a harmonic function, hence the 
name X-harmonic. Moreover, the X-harmonic functions are related to conditioning super-Brownian 
motion in the same way as harmonic functions are related to conditioning Brownian motion; they 
give us the explicit Girsanov transformation to switch from the unconditioned probability law to 
the conditioned probability law. An i?-transform is obtained from P^ by setting 

P?<X) = -^— PlJL {H{X D )Y) 

for Y non-negative and J-" C £>-measurable, where J-qd is the cr-algebra generated by Xf)>,D' C D. 
In his book, Dynkin (2002) suggests a new direction for investigating the solutions of the p.d.e. 

= 2u 2 , namely to explore X-harmonic functions (thinking them as the analogue of harmonic 
functions) and ultimately, to build a Martin boundary theory for this non-linear p.d.e. In this 
case, the Martin boundary is defined as the set of extreme elements of the convex set of all X- 
harmonic functions. Dynkin points out that a concrete understanding of extreme X-harmonic 
functions might yield further insights about the solutions of the p.d.e ^Au = 2u 2 . Since then 
major progress has been made on the study of the solutions of this p.d.e. using other approaches. 
For example Mselati (2004) classified the solutions as cr-moderate in the case of a smooth domain 
(solving a conjecture of Dynkin and Kuznetsov). However, the relationship between X-harmonic 
functions and solutions remains largely unexplored. Dynkin, in a series of papers, has taken concrete 
steps to better formulate and understand extreme X-harmonic functions. Dynkin (2006b) obtains 
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the extreme A-harmonic functions by a limiting procedure from the Radon-Nykodym densities 
H U») = of P »,Xo(dv) = P»(X D G du) with respect to P c ,x D (du) = P C (X D G du). 

Dynkin (2006a) derives a formula for H D (fj.) using diagram description of moments. The functions 
H D will be central to our analysis as well; we will run into them while studying conditionings of 
SBM, and use them to derive results about the structure of conditioned SBM. 

Our goal in this paper is to explore various ways of conditioning a super-Brownian motion. We 
are motivated by the rich structure of the underlying probabilistic objects as well as its potential 
connection to Dynkin's research program on Martin boundary theory of SBM. Here is a summary 
of our contributions: We develop a way of conditioning a super-Brownian motion, which we call 
"conditioning on boundary statistics". The random variables which we condition on are defined 
on an auxiliary probability space, and generated by sampling from the exit measure Xjj. We 
find representations of these conditionings as .ff-transforms of unconditioned SBM. In general, we 
identify the functions H as "extended" X-harmonic. The term "extended" is used because, even 
though we show that these functions satisfy the mean value property (1), in general, we do not know 
whether they are finite for all fi. An example of a boundary statistic is a Poisson random measure 
with intensity /3Xd, where j3 > 0. It turns out that for this kind of conditioning the resulting H- 
transform is through an X-harmonic family of functions studied earlier in the literature. Another 
important and more complex example we studied is conditioning SBM on its exit measure Xd, 
which is what most of our paper is devoted to. We find that the corresponding ^-transform uses 
the family H D {fx) = d p'* D (u), densities first introduced in Dynkin (2006b). This paper shows that 

one can choose a version of this family such that for each u, H u (/i) will satisfy the mean value 
property (1) for each fi. Although this family is claimed to be X-harmonic in Dynkin (2006b), we 
are not aware of any results that actually shows that H v will be a finite function of [i. In this paper 
we shall classify H D as an extended X-harmonic function and leave the question of finiteness to be 
resolved in a different paper, as this will require us to develop analytical bounds on the densities 
of moment measures of SBM. In this paper we will not go beyond describing the probabilistic 
structure of conditioned SBM on its exit measure, as such we are not going to lose much generality 
by stating and proving our results without assuming that H v is finite. Indeed, for our purposes it 
will suffice that H^ D (fi) is finite P^ almost surely, which is true. The family H d (ijl) is of special 
interest because it can be considered as the analogue of the Poisson kernel of D. Also, Dynkin 
(2006b) showed that if H is an extreme X-harmonic function in E then for every fi, and for every 
sequence exhausting E 

H(ji)= limff^Gu) 

Pj^ almost surely. 

The heart of the paper is Theorem 8, which gives a new formula for H D . From this we deduce 
an infinite particle fragmentation system description of P^ = P^" , the conditional law of X in D 
given Xd = v (Thm 11). This is carried out in terms of a particle system, called the backbone 
in Salisbury and Verzani (1999), along which a mass is created uniformly. In the backbone, each 
particle is assigned a measure u at its birth. The spatial motion of the particle is an fa-transform 
of Brownian motion, where h is a potential that depends on u. The measure u represents the 
particle's contribution to the exit measure. At the particle's death two new particles are born 
and u is passed to the newborns by fragmentation into two bits. Here, we use the techniques of 
Salisbury and Verzani (1999) applied to a more general setting. This description connects the theory 
of conditioned super-processes to the growing literature on infinite fragmentation and coalescent 
processes, (see e.g. Bertoin (2006) for a comprehensive exposition). 
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2. Preliminaries. 



2.1. Super- Brownian motion. We will follow Dynkin's definition of super-Brownian motion 
(SBM). Let £ be a domain of M. d and let Me be the positive finite measures on E. A super- 
Brownian motion, [Xd,P^), is a family of random measures (exit measures) and their associated 
probability laws where D is an open subset of a given domain E in W 1 and fi is a finite measure on 
£J with the following properties: 

(a) Exit property: P im (Xd(D) = 0) = 1 for every fj,, and if fi(D) = then P^Xp = /i) = 1. 

(b) Markov property: If Y > is measurable with respect to the cr-algebra J~ c d generated by 
Xrj',D' C D and Z > is measurable with respect to the cr-algebra T^d generated by Xjjh, 
D" D D then 

P ll (YZ) = P ll (YP XD Z) 

(c) Branching property: For any non-negative Borel /, P^(e~^ XD '^) = e~^' VD ^ where 

v D f(y) = -log P v {e-< x °'fi) 

and P y = Ps y . 

(d) Integral equation for the log-Laplace functional: Vnf solves the integral equation 

u + G D (2u 2 ) = K D f 

where Gd and Kp are respectively Green and Poisson operators for Brownian motion in D. In 
other words, if £j is a Brownian motion starting from x, under a probability measure Tl x , then 
KdKz) = ^-xf(Cr D ), where td is the exit time from D. Likewise, GdI(x) = U x (Jq D f(S,t)dt). 

Under certain regularity conditions on D and / (see e.g. Dynkin (2002)), the integral equation 
in (d) is equivalent to the boundary value problem 

-Au = 2u 2 
2 

u(x) = f(x), x G dD. 

Xo represents the exit measure from D, and the first property simply means that Xjj is concen- 
trated on D c , so that exiting is instantaneous if we start outside D. The third property means that 
distinct clumps of initial mass evolve independently. It follows from continuity of Brownian motion 
that X D is supported on dD, if the initial measure is supported on D (see 2. 2. A of Dynkin (2002)). 
The fourth property restricts attention to finite variance branching, and normalizes the branching 
rate. The normalizing factor 2 in front of u 2 is chosen to be consistent with Le Gall (1999) and 
Salisbury and Verzani (1999). 

2.2. Infinite divisibility and Poisson representation. It is well known that Xd has an infinitely 
divisible distribution for each D. This property leads to the construction of a new measure, N x , 
called the super-Brownian excursion law starting from x. Under N x , X evolves as a super-Brownian 
motion, but f$ x will be cr-finite, not a probability. Thus, it is technically more complicated than 
P^. But N x is actually a more basic object heuristically, under which the genealogies are simpler, 
because all the mass starts from a single particle located initially at x. 

In fact P^ can be built up from a Poisson random measure with intensity 6(dx) = J ^ x (dx)/-i-(dx). 
More precisely, let 

U(d X ) =J2 S x i 

i 
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be such a Poisson random measure, where the \ % are random measure valued paths. Then X = 
^2 x l = J X^-idx) is a super-Brownian motion with initial state /i. In terms of Xjj, this yields 
the following formula (see Theorem 5.3.4 of Dynkin (2004)): Let F be a non-negative measurable 
function defined on Me- Then 

oo 

(2) P„{F{X D )) = e-^ M ^F(0) + ^ - e ~ n ^ M ^ / n^dui) ■ ■ ■ K^dujFfa + ■■■ + <) 

n=l ^ 

where IZn is the canonical measure of Xjj with respect to P^ and can be derived from by 
7£ M (^4) = (//, N(.)(X_d G A, 7^ 0)). In other words, we obtain Xq as a superposition of a Poisson 
number of more "basic" exit measures, each descended from a single initial individual. These "basic" 
exit measures arise as the atoms of a Poisson random measure whose characteristic measure is 1Z^. 

Note that a special case of the above representation gives us Vof(x) = N x (l — e~( x °'f'). Note 
further that we will in future write this as 

oo 

P^F(X D )) = ^ - e -^(A^) / n^{d Vl ) ■ ■ ■ K^dvjFiyx + ■■■ + <) 

71=0 n ' ^ 

by taking the convention that the n = term is the first expression in (2). 

The measure N x was first considered by J.F. Le Gall (1999), in his random snake formulation of 
super-Brownian motion. As we follow Dynkin's framework to study our problem, we refer the reader 
to Dynkin (2004) for a systematic account of the theory of the measures N x and their applications. 
Note that the latter has a general branching function ip, but for us, this is taken to be i/j(u) = 2u 2 . 

2.3. Moment measures of super Brownian motion. Among the key tools in our analysis are the 
recursive moment formulae of SBM. The moment measures of SBM are the following measures: 
Let (f>, fi, . . . , f n be positive Borel functions and write / = (/i, . . . , f n ). For C C {1, . . . , n}, let 

(3) n c (^f,x) = n x e-^^n ie c(X D Ji) 

(4) Pc(^/,m) = P^ x ^U ieC (X D ,fi). 

Let K l D and G l D be the Poisson and Green's operator for the operator 

C l = - A — I 
2 

where l(x) = 4Vf)(f)(x). In other words, let £t be a diffusion starting from x, with generator C l under 
a probability measure 11^. Then K l D f(x) = H l x f(^ TD ), where tjj is the exit time from D. Likewise, 

G l D f(x) = ILl(£ D m)dt). 

For C = {i} we have the Palm formula 

(5) nc(<f>,f,x) = K l D fi(x), 

and for general C we have the following recursive formula (see e.g. Theorem 5.1.1 of Dynkin (2004), 
or Lemma 2.6 of Salisbury and Verzani (1999)): 



(6) 



nc{(fr,f,-) = g Yl G l D {4n A {(/)J,-)n c \ A {(t>J,-)) 
AcC,A^9,C 
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(7) 



tt(C) 



Here 7r(C) is the set of partitions of C. These formulas will allow us to construct a variety of 
extended X-harmonic functions of polynomial type. 
We will also need extensions of these formulae for 

® x e-(( x ^)+---+( x ° k >^»n ieC (XD,fi) 

and 

Ppe-« XD ^ + -" + & D k>to))lL ieC (X D ,f i ) 

where D{ C = D. The formulas (5), (6) and (7) give us these quantities when k = 1 and x £ D. 
For k > 2 we find them recursively as follows. Let us put I = {Di, D2, ■ ■ ■ , D k }, 4>i = (4>i, ■ ■ ■ , 4>k), 
u J (s) = N x (l-exp-((X Dl) 0i) + • • • + (X Dh ,</> k ))), I 1 = Au 1 , Dj = DiD- • -nD k and lj = !-{£>,-}. 
We define an operator (</>/, f,x) as follows. For card(J) = 1, 



N x e-« x ^n ieC (X D ,fi) x £ D, 
(8) n&faM = { fi(x) x^D,C = {i} 

x £ D, card(C) > 1 



and for card(J) > 1, 
(9) n&faM 



We let 

(10) P r c (<h,<h,---^k,f,(i) = P M e-(^^ 1 >+-+^'^»n ieC (x i3fc ,/i). 

Fix an / and write Uq(x) for (</>/,/, x) and for p I c ((pii fi A*)- The following formulas tell 

us that the values of n^,(x) on Z?/ can be recovered from its values at the boundary of -D/ and the 
values of the functions n-^(x) on Z?/, with A C C and A / 8,C, using the Poisson and Green's 
operators of C 1 ' for the domain Dj. For C = {i} and card(I) > 1, 

(11) r£(x) = <(«c)W-^D ; , 
and for card(C) > 1 and card(Z) > 1, 

(12) n'dx) = \ Yl G^nWcXA^ + KMcWiXCDi- 

AcC,A^9),C 
Then for /i compactly supported in Dj 

(13) phM = E ^ n k) •••<M, nh r )- 

7T(G) 

Since at the boundary of D/, njj(x) is recursively defined in terms of ?i^(x), we have now a 
complete recursive algorithm to compute nf-,(x), starting with I = {D}, and C = {i}, recursively 
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first building tiq{x) for all C keeping / the same, and then increasing the cardinality of / by 1 and 
repeating the same procedure until the desired cardinality of I is achieved. 

We omit the proof of these formulas and refer the reader to the proof of Theorem 5.1.1 of Dynkin 
(2004). The reader will realize that the argument in that proof works also for the functions Uj{x) 
which are defined as follows. Set Xj = (Xd 1 , . . . , Xu k ), so (Xj, 4>j) = Y^j {Xdj , 4>j)- For for / = {D} 
set 

X/ \ / n x (l-exp-((X I ^ I ) + X l (X Dk ,f l ) + --- + X n (X Dk J n ))) xeD 

Ul{X> 1 <t>k{x) + Xlfl{x) + ■ ■ ■ + Xnfn{x) X £D, 

and recursively for card(/) > 1 as 
u x j{x) 

We note that satisfies for x G Dj 



N x (l - exp -({Xj, 7 ) + Xi{X Dk Ji) + ■■■ + X n {X Dk ,f n ))) xeD! 

u) ix) + <pj (x) x<£Dj,j el,j ^k. 



u) + G lI Dl {2{u)f) = Ki i {u)), 

following the formula 2.11 of Dynkin (2002) Chapter 3, which in turn yields the formulas for 
by differentiating with respect to A. 

2.4. Absolute continuity. The moment formulae together with the Markov property and Poisson 
representation yield an important theorem taken in this form from Theorem 5.3.2 of Dynkin (2004). 
See also Proposition 2.18 of Mselati (2004). Let A4 C D be the space of finite measures compactly 
supported in D. 

Theorem 1. Suppose A^T^ D . Then either P^A) = for all fi G M C D or P^(A) > for all 
fj, G M C D . 

3. Extended X-harmonic functions and conditioning. In the rest of the paper we fix 
D (s E. Following Dynkin (2006b), a non-negative function H : M C D — > [0, co) is called X -harmonic 
in D, if for any D' <s D and any finite measure [i G M C D ,, 

(14) P»(H(X D ,)) = H(ji). 

We will call a non- negative H extended X-harmonic if it satisfies (14) but is not necessarily every- 
where finite. 

We are going to touch upon three different kinds of extended X-harmonic functions, which are 
derived from conditioning SBM on its various boundary statistics. These boundary statistics are 

(a) a Poisson random measure with characteristic measure (3Xd; 

(b) a random variable Z drawn from the probability distribution (^XdI) ^ ^ D ^ ^ an d set equal 
to some given A ^ dD if Xd = 0; 

(c) L(Xr>), where L is a linear map from A4qd to a vector space V (for example, L(fi) = fi or 
L{jj) = (/i, 1) means we condition on Xd or on its total mass). 

Let S be any one of the above statistics. Let X be the state space of S. We will assume that E 
is endowed with a countably generated a-algebra S such that (S,5) is a measurable Luzin space 
(see Dynkin (2006b) for a definition.) For example when S = Xd, S is A4qd, the space of finite 
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measures on dD and S is the <r-algebra in Mod generated by the functions f{y) = v{B), where 
B is a Borel subset of dD. Given Xq = v we let Pg denote the conditional distribution of S. For 
example, Pg{f) equals (u,f)/(u, 1} in the second case (provided v ^ 0), and f(L{y)) in the third. 

Pa denotes a probability measure in which X is an SBM started from fi £ M C D , and in which S 
is then drawn (if necessary) by further sampling. In other words, P^ is a probability defined on the 
fj-field G = 7" cD Vcr{5}. When fi = 5 X we set P M = P x . By construction, P^(f(S) \ T cD ) = Pg D {f)- 
In other words, for any J r C £)-measurable Y we have that 

P lt (f(S)Y)=P fi (lf D (f)Y). 

Likewise, we let P^ 5 and P x 5 denote the marginal distribution of 5 under P„ and P x , so P^ s(f) = 

Let J~cD- = cr{Xf)',D' <s D}. What we want is the conditional law of {Xd>,D' d D} given 
S = s, which should therefore be a transition kernel P^ from S, the state space of S, to the 
J^C-D- measurable functions. More precisely, we will have P^(Y) = P^(Y\S), P^ a.s. for all T<^d- 
measurable Y. The following theorem tells us how we can construct this transition kernel (part 
(d)). The first three statements ((a),(b) and (c)) of this theorem are equivalent to Theorem 1.1 
of Dynkin (2006b) in the case S = Xp. For a general S, we follow Dynkin's proof, with some 
modifications. 

Let us fix a point x £ D. 

Theorem 2. There exists a family of non-negative functions {H x : M. C D h-> R+,s € X} with 
the following properties: 

(a) H x (•) : {s,n) 1 y H x (fi) is measurable and strictly positive; 

(b) For all \i, H x \n) : s 1-4 H^.(fi) is a version of 

(c) For all s, H x (-) is extended X -harmonic in D; 

(d) Define a probability P^ on J- c d- by setting 

(15) I»(Y) = -^—P^YH S X {X D ,)) 

for all D' (s D, and J- cD' -measurable Y , (whenever H x (fx) < 00, and otherwise setting P^ to 
an arbitrary probability measure). Then P^ is a version of the conditional law of X given S 
with respect to P^ for all fx. 

Any two families satisfying the above properties will coincide for P x ^-a.e. sGE. 

Proof. Existence of a family {H x ,s € £} with the first two properties follows from Theorem 
A.l. of Dynkin (2006b) and the absolute continuity of the family {P^s,^, & -^d} with respect to 
P x ,s- Let O be a subdomain compactly contained in D. Then 

(16) P»H' x {Xo) = H S M 

for P Xj 5-a.e. s, V/z € Mq. Dynkin (2006b) proves this when S = Xp , and the proofs for the other 
cases are almost identical to his. Next, we want to construct an extended X-harmonic function H s 
for all s G E. To do this, we choose a countable base O n (w.l.o.g. closed under finite unions), and 
probability measures fi n € A4q , and we let 



R(dr,) = Y,2- n P^(X 0n €d V ). 
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Note that (16) implies 

P tl H s x {X 0n ) = H s M 

for R x P X) g-a.e. (/i, s). By Fubini's theorem we deduce that there exists a P^g-null set M s.t. 
(17) P^H s x {X 0n ) = H s x {n) Vn, for R-a.e. fi G M% Vs € M c . 

For s € W c and /U € A^^,, we choose O n containing the support of ^ and define 

H s M = P,H s x {X Qn ). 

H x (fi) > since this is true of H x , but we cannot rule out H x (fi) = oo. We set H x (fi) to some 
arbitrary positive constant for s £ J\f. The definition of H x ([i) is independent of the choice of O n 
since, if Ok D O n then 

P^H x (X 0k ) = P f ,Px 0n H x (X 0k ) 

= P,K(x 0n ). 

The first equality is due to Markov property. The second equality is due to (17) and the fact that 
P^{Xo n € (•)) is absolutely continuous with respect to R, by Theorem 1. 

Clearly, H x (fi) is measurable and by (16), is a version of dp M '^ (s) for each fi € -M- C D . To show 
that i?|(/i) satisfies property (c), we need to show that H x (fj.) is extended X-harmonic for each 
s € Af c . Let \x and O be s.t. fi G TWq and pick O n s.t. O is compactly contained in O n . Then, by 
definition, 

H S M = P^H s x (X 0n )) 

and 

P,H s x (X ) = P,Px {H s x (X 0n )). 

By the Markov property these two are equal. Now the family {H x , s£S} satisfies properties (a),(b) 
and (c). 

Let us define P* as in (15). It remains to prove that P^ is the desired transition kernel. Let 
D' eg D and Y € T cD >- Then 

P M (/(5)P / f(y)) = J f(s)P*(Y)P, tS (ds) 

/(S) P^) / y ( w )^(^(^))^(^)^,5(^) 



/ /(a) J ff-(X J y(w))P Xj s((is)P M (cL;) 



f( w ) y /( s )p Xo , HiS (d S )p M (^) 

= P M (FP Xjy (/(5))) 

= p^(/(5)y). 

Here we are using the definition of P^,s, the definition of P*, Fubini's theorem, the definition of 
H x , the definition of P v ,s, an d the Markov property of X. 
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Uniqueness follows by a similar argument. Suppose {i^|} s6 s and {i7|} sG s be any two families 
with the properties (a),(b),(c) and (d). Then 

H s x {n) = H*(fi) for P x , s -a.e. s, V/x 

because of property (b). With R(dfj>) as before, there is therefore a -P^s-null set M such that 

H x (fj) = Hxit 1 ) f° r -R-a.e. fi and for s g M. 

Let fj, € M. C D and s M. Choose O n such that /i € M. C Q . Then by absolute continuity and the 
property (c) 

H' x {p) = P»H s x (X 0n ) = P»K(X 0n ) = h s M. 

□ 

In the remainder of §3 we consider the three special cases described above. Our goal is to obtain 
relatively explicit formulae for H x in each case. 

3.1. Conditioning on a Poisson random measure with characteristic measure pXp. Let N be a 
Poisson random variable with mean (Xd,/3). Let Z = {Z\, Z%, ■ ■ ■ } be an i.i.d. sequence of random 
variables from X£)/(Xjj,1). Let 

N 
i=l 

Note that conditioned on Xp, Yg is a Poisson random measure with characteristic measure (3Xrj, 
and that the construction makes sense even if Xq = 0, because then both N and Y equal 0. 

Taking S = Yr, Theorem 2 gives an extended X-harmonic function (which we denote Hx to 
make explicit the dependence on (3) for conditioning on Yg = v. Here v is an atomic measure. We 
let P^ u denote the law of the corresponding conditioned process. In principle this is only uniquely 
defined for a.e. u, but we will find an explicit form that is valid more generally. 

It will be convenient to also define variants of these objects. For any positive integer k, take 
Sk = (Zi, . . . , Zk) if N = k, and Sk = A £ 3D otherwise. Set X^{dz\, . . . , dzk) to be the product 
measure Xj^{dz\) x • • • x Xf)(dzk), and let P^s k ' 3e ^ ne distribution of Sk with respect to P^. In 
other words, for k > 1 and / : (dD) k U {A} — > R such that /(A) = 0, we have 



^pj([ fdxhy 

K! V J{dD)k 



So for any (3 > 0, positive integer k, and any fc-tuple z = {z{\ of elements of dD, Theorem 2 
gives us a family of extended X-harmonic functions {H@ )k,z : z = {z\, . . . , Zk) € (dD) k } such that 

dP & 

(18) H ^( f i) = -^(z 1 ,z 2 ,...,z k ). 

dP x,s k 
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Let lp = 4Vb(/3). We have P^YpiO} = P^{e~( XD ^) = and P X) y„{0} = e~ l ^ x \ therefore 

we find H@' v for v = simply by the ratio 

by Theorem 2. 

Let ( > be a bounded Borel function on D. For i £ D, we let m l x (dz) = H l x (^ TD € efe) denote 
harmonic measure on dD for the operator L'. Then m l x and are mutually absolutely continuous, 
for x,y € D. (This is a well known fact; however, for the curious reader, here is a quick argument 
why it is true. Let D' be a smooth domain compactly contained in D, and x,y € D'. Let D / be 
the harmonic measure on dD' . If A is a Borel subset of dD and m^-A) = 0, because of the strong 
Markov property and the fact that m l x D , is equivalent to the surface measure 773/ on dD', we have 
that m l z (A) = for 773/ almost all 2. This implies m^A) = 0, again due to the strong Markov 
property and the fact that m l y D , ~ 7/)'. ) 

Let 

, dm,, 

»«».«) = 3^ w 

denote the density. If -D were sufficiently regular, this would be a version of the Martin kernel for 
the operator L l , but we make no such regularity assumptions at this point. We take k l to be a 
jointly measurable version of this density that is harmonic in y, for each 2 E dD. One can construct 
k l in a similar way as in Theorem 2. That is, we start with a family {&((.), z),z € 3-D} such that k 

is measurable as a function of (y, z), and for fixed y, k(y, (.)) is a version of j^f~( z )- The existence 
of such a family follows from Theorem A.l. of Dynkin (2006b) and the absolute continuity of the 
family {m l y ,y € D} with respect to m x . Then we take a sequence D n <£ D exhausting D, and let 

k l (y,z) = ILy(k(£ TDn , 2)) for y £ D n . Then we prove that k ((.),z) is well defined, and harmonic 
for all z except on an m^-null set A/", on which we set k l to be an arbitrary constant. We omit the 
details as the arguments are very similar to those in the proof of Theorem 2. 
In the case I = we write m x (dz) = m x (dz) and k x (y, z) = k x (y, z). 

The particular case of interest is I = lp = 4Vd(/3). Suppose that k > 1 and that z±, . . . , z\. € dD. 
For C C K = {1, . . . ,k} recursively define 



Finally, set 



forC = {0 

k \YjAcc$^a^c g d(^PaPc\a)' for 1^1 > l - 



01,... ,z k ) 



where the sum ranges over all partitions {C±, . . . , C r } oi K. 

In the following theorem we use the convention that H^ ,0,z = H^'°. 

Theorem 3. Let D <s E, and (3 > and x G D. Then 

(a) Hx = Hx ,k ' z for P x ,Ys' almost all finite atomic measures v, where k and z are such that 

k 

(20) u(dx) = Y / S Zi (dx). 

1 
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(b) For (irix) k -a.e. (z±,..., Zf.), for all /i € M C D 



(2i) - 



0,k,Z(,,\ — \ Zl '---l Z k) 
Px {Zi,...,Z k ) 



(c) If D is smooth, then in fact p^,' (z±, . . . , Zk) < oo for all /i G M C D whenever zi, . . . , z^ are 
distinct. 

Proof, (a) g k (f) remains unchanged if we permute the arguments of /. Thus we can choose 
the densities Hx' k,z (p) to be both X-harmonic and invariant under permutations of the 2j. A 
simple way to confirm this is to replace an X-harmonic choice of Hx ,k,z (p) by ^ Hx' k ' a (ft) , 
where the sum is over permutations a. The latter is still X-harmonic, and a version of the density 
S k /dPx ^ s a ^ so clearly invariant under permutations. 

For a finite atomic measure u, all of whose atoms have mass 1, find k and z±, . . . , z k such that 
(20) holds. Then define 

Note that Hx{p) is well defined, since Hx ,k ' z depends only on z k := (z\,...,Zk) and is invari- 
ant under permuting z k . (Note, if two sequences z and z satisfy (20), then z k and z k must be 
permutations of each other.) 

Let / be a function defined on the space of finite atomic measures. If v = 2?=i write fk{z) 
for f(u) . To finish the proof it is enough to observe 

P X ,Y P (HM (/*)/(•)) = P x (H^(n)f(Yp)) 



k=l 



h£**MMz) XDidz) 



(x D ,i) k 

WOi/W + E^M 6 "^ / H^' k ' z ^)fk(z)X k D (dz) 
k=i ' J 

oo 

p»xAo}f(o) + Y, p U( H % M ' ] (riM-)) 

k=l 

oo 



k=l 



^{o}/(o) + E^(e"^^^#i { x^ 0} / AM^$ 



(b) Define Hx' k ' z (p) to be the right hand side of (21). Following an argument of Dynkin (2004), 
chapter 5 one can show that is the density of P^ $ k with respect to {m x 3 ) k . The argument uses 
the moment formulae (5), (6), (7) and then pulls k factors of harmonic measure out of the resulting 
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expressions, leaving the densities k l £ behind. It follows that Hx ,k,z {p) is a version of the Radon 
Nikodym derivative in (18). The finiteness condition for p^' k follows immediately. 

Furthermore, Hx ,k ' z falls in the family of functions considered in Salisbury and Verzani (1999), 
Theorem 3.1 therefore shows that Hx' k,z is X-harmonic. Thus Hx' k,z = Hx' k,z for (m/) fc -a.e. z, 
which is the sense up to which Hx' z is well defined. 

(c) The argument for (c) is a straightforward modification of the estimates used in Theorem 5.3 of 
Salisbury and Verzani (1999). 

□ 

Remarks: (i). The conclusion is that we have obtained an explicit formula for Hx ,u (pi). The abstract 
definition of this X-harmonic function was valid only up to an unspecified null set of z/'s, whereas 
the canonical expression we have obtained is well defined as long as v is a finite atomic measure, 
all of whose atoms have mass 1 (assuming that D is smooth). 

(ii) . The arguments of this section would work equally well for conditioning on the value of a 
Poisson random measure with characteristic measure /3(x)X£)(dx), where (3(x) is now a bounded 
measurable function on dD. 

(iii) . If D is smooth, then instead of taking k l x (y, z) to be the density of m l y with respect to m x , we 
could use the Poisson kernel in its place, and get a similar result. In other words, we could take the 
density of m l y with respect to the surface measure 7 on dD, rather than the density with respect 
to m x . 

(iv) H^ ,k ' z falls in the family of X-harmonic functions considered in Salisbury and Verzani (1999). 
This family of X-harmonic functions are characterized by a function g, and /3 4fir -harmonic functions 

In our example the function g is up = Vd(3, and the harmonic functions Uj are the 
functions k l x(-,Zi). In Salisbury and Verzani (1999) it is shown that for D Lipschitz of dimension 
d > 4, g = 0, and Vi = k x ( m , £j) where z\, . . . ,z^ are distinct points chosen on the boundary, the 
resulting X-harmonic function corresponds to conditioning SBM to hit the points Zj. The same 
argument would work in dimension d = 3, at least when D is smooth. 

3.2. Conditioning on a r.v. Z sampled from measure ^x^i) • R ecau that the random variable Z 

is drawn from the probability distribution /Xd,i) ^ ^ D ^ ^> an d set equal to some given A ^ dD 
if Xjj = 0. Applying Theorem 2 gives us a family of extended X-harmonic functions 

x dp x y ' 

indexed by points z of {A} U dD. We denote the law of the corresponding conditional process by 

pz 

Recall that £j is a Brownian motion under H y . For z G dD, we let H z be a probability under 
which £ t is a k x (-, z)-transform of Brownian motion. (Recall k x ((-), z) = k x ((-), z)). In other words, 

n*(/(&),* < r D ) = —L_u y (f^ t )k x (^z),t < TD ) 

KxyUi z ) 

for every bounded measurable /. 

The following result establishes a concrete formula for H z that is defined for m^-a.e. z G dD 
when D is a general domain. When D is smooth, the same argument as in the previous section 
gives a canonical version, defined for all z G dD. 
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Theorem 4. Let D <s E and x e D. Then for m x (dz)- almost all z G dD, H^(fi) < oo, and 

= ^^U-,z)W { . ) {e-^)))e~^)d^ 
( } x(P) ~ | oo n|( e -^))e-^Wd/3 

for every \x, where up = VoP and the random variable (p(up) is defined as 
(23) <j,(up) = 4 / Uf) (Zt)dt. 



Proof. We first find the Radon Nikodym derivative of PuZ w.r.t. the harmonic measure m x {dz) 
on the boundary of D. We observe that 

(24) p fM f) - ^(f^W»>) = - f 

Note that the above derivative equals when Xd = 0. By the branching property, 

(25) p ( e -Mx D ,f)-p(x D ,i)j _ e -(jj,,ux}+p) 

where 

uxf+p = V D (Xf + P) = N (0 (l - e -^M+P)y 
Taking the derivative of the right side of (25), and evaluating at A = we get 

(26) ^(f^ 1 *****) = / 00 <^ N (-)«^»/> e " /,(XD,1> )> e "^ )d / 3 - 

Differentiation under the integral sign is easily justified. By the Palm formula, 

N y ((X D ,f)e-^ x °V) = n y (f(£ TD )e-<f>^) = [ Ul(e-^)f(z)m y (dz) 

JdD 

= [ U z y (e-^)k x (y,z)f(z)m x (dz). 
JdD 

So, 



Hence 



( Mj N ( . ) ((X D ,/)e-^^- 1 ))) = / (^k x (;z)Ul ) (e~^))f(z)m x (dz). 

JdD 

(/x,N ( . ) ((X I? ,/)e-^^' 1 )))e-^>d/3 

/(*)( ^ \ t i,k x (-,z)Ul ) (e-^))e-^d^ m x (dz). 



1 3D v JO 

Therefore, both P^z and P x> z have densities with respect to m x (dz), given by 



and 

/"OO 

n£( e -*(«/j)) e -«/>(*)d0 
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respectively. The ratio of these two is a version of the desired Radon Nikodym derivative. To show 
that it equals for almost all z, it simply remains to show that it is extended X-harmonic. 

The denominator is simply a normalizing factor, so consider the numerator. For lp = &up, it is 
known (see Theorem 1.1 of Salisbury and Verzani (1999)) that pL h4 (fi,v)e~ ( ~ fM ' u ^ is X-harmonic 
whenever v is /^-harmonic. And in our case, we have that (fi, k x {-, z)W^{e-^ u ^))e-^ u ^ = 

(/i, k l x (-, z))e~^ t ' u ^ as required. 

Now we show that H*((a) < oo for all /x, for m^-almost all z. Let hq € M. be fixed. Since H*(no) 
is a density, there exist a m x null set -B C dD s.t. for z & B, H*(fJ>o) is finite all z € -B c . Let 
/x € Ai c D and choose D' <s D such that both /i and hq are compactly supported in D', and assume 
that D' is smooth. Let p M be the measure defined on 3D' by Pfl (f) = P^e-^o'^p) (X D >, /)). Then 
our analysis in Section 3.1 gives us that and p^ are equivalent and the Radon Nikodym density 
is 

dp» , v = (/x,fc^((-),y)) 
d Pw [y> (f, ,~kh((.),y)) 

where k P(u,y) is the Poisson kernel of D' for the operator C l P . Note that this density is bounded 
by a constant C(/i,/io) since k lf3 {u,y) is harmonic in u on the support of \x and Hence 

< C( Ai ,/ i0 ) J P M() (X D ,,^(-,z))e-^" u ^ 
= C(/i, Mo )(Mo,^(-,z))e-^). 

It follows that H*([/,) < C(/i, (1q)H*([1q) < oo for all z € B, hence the proof is complete. 

□ 

3.3. Conditioning on a linear function of Xe>. Let L be a linear and measurable map from the 
linear cone of positive finite measures Mod on dD to a Luzin measurable space (V, V) where V is a 
vector space and V is countably generated. Let V + be the image of Mod-, and write V* = V+\ {0}. 
Assume that L/i = implies /i = 0. 

Let T n be the map V+ ->■ V+ defined by: 

T n («i, ... ,U n ) !->•(«! + « 2 H h « n ,«l, • ■ • ,«n-l), 

and for ,4 € V and icF, let 

N xM x D )(A) = n x (L(X D ) G A,X D ^ 0). 
We fix x 6 D as before and define a reference measure R x on y* by 

j^(i) = P Iii(Xfl) (i,I^0). 

Its total mass is r Xj o = P x ,l{x d ){Y*) = 1 — e - u ( x ) ) where u{x) = —\ogP x (X£> = 0). Note that 
u = lim^^oo Vd((3). Throughout this section we will assume that 

(27) D is a bounded domain, all of whose boundary points are regular. 
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This holds, for example, if the boundary of D is smooth. Under assumption (27), Vd{@) is the 
unique solution of 

-Ait = 2u 2 
2 

on D with u = (3 on dD. (Proposition 8. 2.1.B of Dynkin (2002)). Because of this and the comparison 
principle (Proposition 8. 2.1.H of Dynkin (2002)), Vd(/3) is also the maximal solution of \Au = 2u 2 
on D bounded by j3. We will need this in the proof of Lemma 5. 

By Theorem 2 we know the existence of a family of extended X-harmonic functions {H£,v G V} 
such that H%(/j.) = dP^ux D )/dP x ux D ){ v )- I n this section we are going to find a more explicit 
formula for this family. 

Lemma 5. Assume (27). There exists a family of functions {"y x ,v '■ D t- > (0,oo),v G V*} such 
that the mapping (v,y) i-> J x ,v(y) * s measurable and for all y G D 

(28) N yMXD) {dv)= lx , v {y)R x {dv). 

In addition, there exists a measurable kernel K XiTl (v; dv\, dv 2 -, ■ ■ ■ ,dv n —i) from V* to (U*) n_1 , such 
that 

(29) R% o T~ l (dv,dv 1 ,dv 2 , ■ ■ ■ ,dv n -i) = K Xtn (v;dvi,dv 2 , ■ ■ ■ ,dv n -\)R(dv). 
Moreover K x n (v, ■) is a strictly positive measure, for R-a.e. v. 

Proof. Recall that 

P*HX D )(A) = P^HXd) e A). 

Because all P^Xd G •) are equivalent so are the P ll ux D )i as are their restriction to V* . Thus all 
Pfi,L(x D ) (when restricted to V*) are equivalent to R x . Moreover, since 

N yM x D M) = N y (L(X D ) £A,X D ^ 0)) = N y (Px D ,(L(X D ) G A, X D ^ 0)) 

for all y G D and D' <e D such that y G D' , N y ^x D ) is a ^ so equivalent to R x for all y G D. By 
Theorem A.l. of Dynkin (2006b) we get a family of functions {j x ,v : Z? 1 — >- [0, 00), v G V*} such that 
the mapping (v,y) i-> J x>v (y) is measurable and (28) holds. Clearly such a ^ x ,v(y) can be chosen 
strictly positive since N y ux D ) anc ^ ^ x are equivalent. 

It will be convenient for the proof to write R x {dv) = r x $R x {dv), where R x is a probability 
measure. In other words, R x (dv) = P x (L(Xd) G dv \ Xd ^ 0). Note that with this choice of R x , 
R™ o T~ 1 (dv, dvi, dv 2 , ■ ■ ■ , dv n -i) is the joint distribution of (Vi + • • • + V n , V\, . . . , V n -\) where Vj 
are independent random variables with distribution R x . Let R x ^ n be the marginal distribution of 
V\ + • • • + V n , where the Vi are as above. The following decomposition is then immediate: 

i?" o T~ x (dv, dv 1 ,d,V2, • • • , dv n -i) = K x>n {v\ dv 1 ,dv 2 , dv n -i)R X)n {dv) 

where K x>n (v; dvi,dv 2 , ■ ■ ■ , dv n -±) is the conditional probability kernel for (Vi, . . . , V n -i) given V\ + 
VV n .' 

We now show that R XjU is absolutely continuous with respect to R x . Let Xp, . . . ,Xjj be n 
independent realizations of the exit measure under the law P x . Then the distribution of Xj^ + • • • + 
Xp is given by the P n s x distribution of Xd. 
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Let F be s.t. R X (F) = 0, i.e. 

(30) P x (F(L(X d ))1 {Xd ^ 0} ) = 0. 
Because P x and P n g x are absolutely continuous, (30) implies 

(31) P n5x (F(L(X D ))l {XDm ) = 0. 
Since 

P nSx (F(L(X D ))l {XDm ) = (P x ) n (F(L(Xh + ---+XZ))l {xh+ ,,, +X n^ } ) 

> {P x ) n {F{L{X l D ) + ■■■ + L(X n D ))l {xhm ■ ■ • 1 {X ^ 0} ) 
= rl R x , n (F), 



this implies that R Xtn (F) = 0, so indeed, R X)n is absolutely continuous with respect to R. 
If h x (v) is the Radon Nikodym derivative of R x ^ n with respect to R x , we get (29) with 

Kx,n(v;dv 1 ,dv 2 , ■ ■ ■ ,d« n -i) = K x ^(v;dv 1 ,dv 2 , . . . ,du n _i)/i^(u). 

It remains only to show that K xn (v, •) is strictly positive. Because K x ,n( v , •) is ; this amounts 
to showing the converse to the absolute continuity result above, namely that R x is absolutely 
continuous with respect to R x ,n- 

Our approach is to use the Poisson representation, as in the absolute continuity argument in 
Dynkin (2004). Suppose R XjTl (F) = and < F < 1. The Poisson representation gives that 

P^F(L(X D )),X D ^0) = Y,^TT- / / W + - + %R(^e^i,lD^0) 

k=l ' J J 



k=l 

where 



• • • ^x n ( x D G dv k ,X D ^ 0)fi(dx 1 ) ■ ■ ■ n(dx k ) 
2^ — ^ — / fk(xi, ■ ■ ■ ,x k )ii(dxi) ■ ■ -n(dx k ), 



f k ( Xl , ...,x k ) = J F(L(u!) + ■■■ + L(v k ))N Xl (X D e dut,X D + 0) • • ■ N Xk (X D e di/ fe ,X D ^ 0). 



Let D m d L> such that x € D m and -D m f -D. Then 

fUF(L(X D )),X D ^ 0) = P x (P Xl , m (F(L(X D )),XD # 0)) 

00 -<x Dm ,«) 



P *(z_> fcj / fk(xi, ■ ■ ■ ,x k )X Dm (dx 1 ) ■ ■ ■ X Dm (dx k )). 

k=i ' J 



There is a similar Poisson representation for R Xtn (F), involving a sum of integrals of the f k for 
k > n. Since R x ^ n (F) = 0, we conclude that for each k > n there are xi, . .. ,x k S D such that 
f k (x±, . . . , x k ) = 0. By absolute continuity, we conclude that f k (x±, . . . , x k ) = for every x\, . . . , x^. 
Since F < 1 we obtain the bound 

fc 

/(xi, . . . ,x k ) < Y[N Xj (X D / 0) = u(xi) ■ ■ ■ u(x k ). 

3=1 
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Therefore 

P X (F(L(X D )),X D + 0) < Y^P*( jrj (X Dm ,u) k ). 

k=l 

The result will follow once we argue that all terms e'^ ^'^ (X Dm ,u) k converge to P x -a.s. as 
we let m — > oo, by the dominated convergence theorem since these terms are bounded. To show 
that e-( XD ™^(X Dm ,u) k -> 0, it is enough to show that the stochastic boundary value of u (that 
is, lim m _ s , 0O (Xo m , u)) is or oo, i-^-a.s. Let Z$ be the stochastic boundary value of the constant 
function j3. The sequence (Xjj m ,/3) is a uniformly integrable martingale with respect to P y for all 
y £ D. That it is a martingale follows because the constant function /3 is harmonic. The uniform 
integrability follows because this martingale is square integrable. Indeed, by the moment formula 
(7), P y ((X Dm ,p)) 2 = Gz5 m (4/3 2 )(y) + /3 < G D (4/3 2 )(y) +/3 < oo, since D is bounded. It follows that 
the log-potential of Z is the maximal solution of = 2v? on D bounded by /?, which is VdP 
as we argued at the beginning of this section. (See sections 9.2.1 and 9.2.2 of Dynkin (2002)). So 
V D /3(x) = - log P x e~ z ». Also note that = f3Zi, and Zp t Z where 



(32) Z 



if Z x = 
oo if Z x > 



By the dominated convergence theorem P x e~ z = \\m.^^ r00 P x e~ z ^ = lini3_j. 00 e~ v °^^ = e~ u ^ x \ 
so u is the log potential of Z, therefore Z is the stochastic boundary value of u. Since Z is or oo 
P x -a.s., the proof is complete. 

□ 

Lemma 6. Assume (27). Let (K Xtn ) n >2 be a sequence of transition kernels satisfying (29) for 
n > 2. Define 

(33) K X)n (v;dvi, . . .,dv n ) := K n {v;dvi, . . .,dv n -i) x 

"v— (uiH h^ji-l 

)(dv n ). 

Then for R x -almost all v £ V* , the following holds for all n > 2, 1 < r < n, and any partition 
Ci, . . . , C r of {1, ... , n}, where rii = 

K x ,n(v,dvi, . . . ,dv n ) = / K x>r (v;dvi, . . . ,dv r ) K x ,m(vi; dv Ci ) ■ 

^ i=l 

Proof. Since V is countably generated, so is V n (the product <r-field), and therefore for each 
n > 1, there exists a sequence of nonnegative Borel measurable functions generating V n . 

It will suffice to show for any k, j > 1, n > 2, 1 < r < n, and any given partition Ci, . . . ,C r of 
{1, . . . , n} that 



/fe (vi, . . . jtvO-K^nCMt'i, . . • ,dv n ) 



Rx(dv) 



f(v) 



K x , r (v,dvi,...,dv r ) / n r i=1 K Xtni (vi,dv Cl )fk(vi, . . . ,v n ) 



R(dv) 



Let 



y n-^^n^Ui, dv Ci )fk {Vl, ■ ■ ■ , V r , 
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Then, 



fj(v) J K XjT (v, dvi,..., dv r ) J U r i=1 K Xin .(vi, dv c Jf£{vi, . . . ,v n ) 

= J f}(vi-\ hv r )f%(v 1 ,...,v r )R x (dv 1 ,...,dv r ) 

= J fjih + ■■■ + v r )fk(vi, ...,v n ) [Yl r i=l K x ^ i {vi,dvc l )R x {dv l )\ 

= J /j(Sd" m H \-J: Cr Vm)fk(vi, ■ ■ ■ ,v n )U r i=l U c% Rx(dv m ) 

fj{vi H h v n )f% . . . , v n )R x (dvi, dv n ) 



R{dv) 



fj(v) 



/£ fa, . . . , v n )K xn (v, dvi,..., dv n ) 



Rx(dv). 



□ 



Lemma 7. Assume (27). For any B C V* with R x {Bq) = 0, there exists B C V* with 
R X (B C ) = 0, such that B C Bo and for all n > 2 and v G B, we will have (v±, . . . ,v n ) G B n , 
K x>n {v,dv\, . . .,dv n ) a.s. 

Proof. Define recursively B m ,m > 1 as 

oo „ n 

B m = {v G B m _i 1 B^ n _ 1 (vi)K x>n (v, dv\, . . . , dv n ) = 0} . 

n=2^ i=l 

Then i? x (f?£J = 0. Because, R x (Bq) = 0, and assuming R x (B!f n _ 1 ) = 0, we have that 

» oo « n oo « n 

1^2 ^lBz n _ 1 (vi)K x ,n(v,dvi,...,dv n )R x (dv) = ^ / ^ (v^R^dvi, . . . , cfc n ) 

n=2^ i=l n=2^ i=l 

oo 



n=2 



which implies R x {v : Y^=2 I Ya=i _ { v i)^x,n{v, dv%, . . . , di> n ) ^ 0} = 0. This implies that 
= since B c m C {v : f Eti 1 ^ t (vi)K x , n (v, dv u . . . , dv n ) + 0}. Let B = D% =0 B m . 
Clearly R X (B C ) = 0andBcB . 

Now, if u G B then for any m, ^^L 2 / SiLi 1 B^ (vi)K x ^ n (v, dv\ , . . . , (fo n ) = 0. By the monotone 
convergence theorem this implies that Y^=2 I Si=i ^-B c (vi)K x n (v, dvi, . . . , dv n ) = 0. Hence for all 
n, K Xjn (v, dv\, . . . , dv n ) almost all (v±, . . . , u n ) is in B. □ 

So far Tir,'!; (y) is any jointly measurable version of the Radon Nikodym derivative of N yt x D with 
respect to R x . In the following theorem, we refine this choice and find a formula for the family of 
extended X harmonic functions {H x ,v G V*} corresponding to dP fl x,(x D )/dP Xt L(x D )( v )- 

Theorem 8. Assume (27). There exists a version of {j XtV ,v G V*} of Lemma 5 such that for 
R x almost every v G V* 



(34) 



j x , v (y) = N y (H x '(X D ,)) 
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for every y, and D' d D such that y 6 D' , where 

' e ~{^u)+u(x) ^ if v — 

(35) H%{n) = I g-^^^ + ^j^e-^)^^;^,...,^^))^ 

X(Mj7«i> ••■(/*, 7*,tt n _i>(At>7!r,t»-(tii+-+«n-l)>» { ! v + ^ 

Moreover, with such a version of ^ x ,v, {H%,v £ V*} defined by (35) is extended X-harmonic for 
R x -almost all v. For fixed fi, H x (fJ,) is a version of dP^ t L(x D )/dP Xt L(x D )( v )- 

PROOF. Let {j x ,v, v € V*} be a family of functions satisfying the properties of Lemma 5. Define 
H% by formula (35) with this j XtV in place of Jx,v Since j XtV > 0, it follows that H x > 0. We will 
first show that for each H x ([i) is a version of dP^uxn) I dP x ,L(x D )( v ) ■ Let F be a non-negative 
Borel function on V+. 

P xM X D )(dv)F(v)H v M 

= P xMXD) ({0})mS°M + j F(v)Hl(»)R x {dv) 

= e - u{x) F(Q)e-^ u)+ui - x) +e~ M J F(v)(fi,% v )R x (dv) + 
~ e -(Mi«> r r 

+ 2^ , / / <Mi7*,Wi) (A*>7x,» ft -l)<A*>7!r,t;-(wi+».+«n_i)> X 



n=2 



xF^iC^^;*!,. . . ,dv n -i)R x {dv) 

-(fi,u) 



e~ M F{0) + — / F(vi + ■■■ + v n )(fi, % m ) ■■■{(!, %, Vn )R x {d Vl ) ■ ■ ■ R x {dv n ) 

n=l n ' ^ 

e -M F (0) + £ ^— 1- / F(vi + ■■■ + v n )((i, N { . )MXD) (d Vl )) ■ - - </i, iV ( . ))i(XD )(^ n )) 



n! 

rc=l 



oo —(fj,,u) P 

\ F ( L ( v i) + ■■■ + L(u n ))n^du!) ■ ■ ■ K^dvn) 



n=0 

= P, L (F(L(X D )) = Fi , MXd) (F) 

where we are using the Poisson representation. We choose countable base O n (w.l.o.g. closed under 
finite unions) and consider the measure R on M C D defined as in the proof of Theorem 2. The 
argument in that proof tells us that there exists a P^pf^-null set Af s.t. 

(36) P^H v x {X 0n ) = H l x {fi) Vn, for i?-a.e. fx € M C D , \/v € Af c . 

Without loss of generality we may assume that for v G Af c , the statement of Lemma 6 holds, and 
moreover by Lemma 7, K n (v, dvi, . . . , dv n ) almost surely (v±, . . . , v n ) € (Af c ) n for all n. For v € Af c 
and y £ D, we choose O n containing the support of y and define 

lx , v {y)=n y {H v x {X 0n )). 

Since H x > it follows that j XiV (y) € (0, oo] for every y. We set -f XiV to some arbitrary constant 
> for v € TV. The definition of ^y X)V is independent of the choice of O n since, if Ok D O n then 

N y HZ(X 0k ) = N y Px 0n S s x (X 0k ) 

= N y H s x (X 0n ). 
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The first equality is due to Markov property. The second equality is due to (36) and the fact that 
N y (Xo n € (•)) is absolutely continuous with respect to R. 



We now show that for fixed y, this gives a version of — V dR XD) • ^ V ^ ® n then 

(v)N yMXo) {dv) = n y {F(L(X D ))l {XDm ) 

= N y {P Xon (F(L(X D ))l {XDm )) 



N„(/ Hl{X 0n )F{v)R x {dv)) 
F(v)R x (dv)N y (H^X 0n )) 



F(v)j x , v (y)R x (dv). 



Let H x be defined by formula (35). Then we know that for fixed fi, H x {ii) is a version oldP^^ XD ^/dP x ^ XD ^{v) 
Now we are going to show that for v € M c , H x (fj.) = P fl (H x (Xo k ))- To simplify the notation, we 
drop the basepoint x from K x ^ n , j XiV and j X)V for the remainder of the proof. 



P^Hl(X 0k 



~ r P -(Xo h ,u) 



n=l 



71 



-Ui =1 {X 0k ,j Vi }K n (v;dv!,..., dv n ) 



oo „ 

J P ^{ e ~ {X ° k ' U) K=i(Xo k ,7v i ))K n (v-,dv 1 ,...,dv n ) 



n=l 

oo 



n=l n ' •* 



Ufa, N (0 (er<*>k > u >11 jeGi (Xo k , % ) 

7r(ri) i=l 



K n (v;dvi,. . .,dv n ) 



oo r 



e -{M,V 0fc («)) 

n=l "" n (n) u i=l 

K r (v; dv\, . . . , v r ) 



1 n f 

S^iS / ^r(«;«l,-..,«r) 

1 " 1 ^ 



(■) 



n=l r=l 

-<X 0fc ,«> 



ml ■ ■ ■ n r \r\ 



i=i 



-<M,Vb fc («» 



-(M,Vo,(«» 



-(Xo fc ,u> 



n i6Ci (X 0fe , 7 Uj )lf nj («i ; dv Ci ) 



U jeCi {X 0k , Ivj )K rH (vi ; dv Ci ) 
£ / - ] K r (v-dv u ...,dv r )]\ ^T^' N (-) ( e 

r=l ' i=lrii=l l ' ^ 

oo ,. ^ r | °° 1 /" 

£ / -ir r (^,---,^)II<^ N (-) 2 — e~ <Jf ° fc,u> / ' I^o,,^,)^^;^) 

r=l J r - i=l W«=l 7 ieCi 



r=l ^ r ' i=l 



where the last line follows because by our assumption that if v € M c , then K r (v , dvi, . . . , dv r 
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almost surely (v±, . . . , v r ) G (M c ) n and by construction for each V{ 

oo . 

ivM = n„ —, e ~ {x ° k,u) / n^c^Xo,,^.)^^;^). 

Now we show that 7 X) „ satisfies equation (34). Let y £ D D. Then for some n, L> <s O n . By 
the definition of 7^ and the Markov property, 

N„(flJ(X 5 )) = N J/ (Px 6 ^(^oJ)=N 2/ ^(X J) = 7^(2/) 

as desired. 

Replacing N y with in the above argument we have that H% is extended X-harmonic for all 
v € Af c . It is also clear that H v x is X-harmonic for v = 0. Hence the proof is complete. 

□ 

Remarks: 

(i) Although Theorem 8 gives a workable form for H^(fi), it does not remove all ambiguity in 
the choice of flj, since 7e,d and K x (v'i ■) are only well defined for a.e. i>. Ideally we would like to 
prove continuity properties of these objects in i> as well, that would then specify them uniquely. 
But we have not succeeded in doing that. In subsequent sections we will however be able to clarify 
the structure of these objects, and show how they determine the behaviour of the P£-transformed 
super-Brownian motion. We do have regularity in y. In particular, we will soon see that the version 
of lx,v(y) given by Theorem 8 is already lower semi continuous in y. 

(ii) Note that in sections 3.1 and 3.2 without much difficulty we were able to show that the 
corresponding extended X-harmonic functions are finite, therefore X-harmonic. Finiteness is harder 
to prove for the extended X-harmonic functions of this section because it requires analytical bounds 
on the Radon-Nikodym densities of n-order moment measures of SBM for all n > 1. We hope to 
pursue such bounds in a subsequent paper. 

(iii) An important case is when L(Xd) = Xjj. The corresponding extended X-harmonic function 
H£ can be thought as the analogue of the Martin kernel. 

(iv) A more tractable application should be the case L(Xd) = (Xd, 1), where we condition on 
the total mass. In that case, j x ,v(y) is a function of only finite dimensional variables v 6 [0, 00) and 
y € D. We hope to explore this example further in a subsequent paper. 

(v) An interesting direction is to explore the relationship between Hx U o£ Section 3.1 and - In 
particular, with (3 = n, what happens to H n,Vn as n — > 00, if {vn}n>i is a sequence of finite atomic 
measures such that n~^v n converges to a finite measure v on the boundary? One can show that 

weakly almost surely, in a sense, and use this to investigate whether H% is extreme. We will pursue 
this direction in a subsequent paper. 

4. Fragmentation system description of PJ£. The results of this section apply in general to 
the conditional law given L(Xd) = v. For simplicity, however, we will carry out the computations 
for the case L{y) = v. So in this section V* is the set of positive finite measures on dD. 

As in Section 3.3, we assume (27), ie that D is regular. From Section 3.3, recall that u(x) = 
— \ogP x {Xn = 0) = Vb(oo) and u is a solution of the boundary value problem 



-Au = 2u 2 
2 
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on D and u = oo on dD. We also have 



(37) u(x) = n x (X D > 0) = lim V D (3(x). 

/3->oo 

Let j XjU and H^. be as constructed in 3.3. Let N C V* be the -R x -null set such that for £ AA C , 
7 X)J/ and i?^ satisfy the system of equations (34) and (35), H% is extended X-harmonic, and the 
kernels K x ^ n are strictly positive and satisfy the decomposition property of Lemma 6. By Lemma 7 
we may assume that if v G M c then K x ^ n {u, du\, . . . , df n )-almost all (u±, . . . , u n ) are in (Af c ) n . We 
fix v € A/" c . Let y be a point in D such that ~f X)U {y) < oo. Recall that ^f X) u{y) > by construction. 
If y € D' d -D we may then define a change of measure by 

W y {Z) = -J—N y (ZHZ(X D ,)) 
ix,i/\y) 

for positive J^ D i -measurable Z. (J^pi is defined as the a-algebra generated by {Xo,y € O C I?'}) 
Since -£f£ is extended X-harmonic, is defined consistently on J r ^- D _, and we have that is a 
probability law because of equation (34). 

In the remainder of this section we turn to the problem of giving an explicit probabilistic con- 
struction of in terms of a backbone along which unconditioned mass is created. We do this in 
two steps. Let H% = e~^^ (/x, 7 Z)t ,), and for n > 2 let 

1 7x,U-(ui-\ h^n-l)/' 

Then H% = Xm>i H xn . We will first use the recursive moment formula to establish an inductive 
relationship for H% n , and then compare this to the inductive relationship coming from the first 
branch of the backbone. 

We will make use of a stochastic process £t under various measures 11^ or n^ u . In either case we 
use the shorthand 

A/f'(0) = e -/o 4 ^ s (l-e-^"«)d, 

where D' <s= D and let tq/ be the exit time of £ from D'. We similarly let 



Mi 



where / = {D\, . . . , D^} is such that each Dj (s D, cpj = (fa, . . . , 0*.), Xj = (Xm, ■ ■ ■ , ^D k ), so 
(X T , fa) = "£AX D , fa). Also write D I = D 1 n---HD l 



fc- 



Ill the rest of the paper we will denote K x ^[y^ dv') by K x [y, dv'). Let / = {Di, . . . , D^} where 
C D k = D' m L 
For card(I) = 1, 



Di C -Dfc = D' (s D. Define a family of operators Ny' u,n as follows: 



n y (e^ x D'^H» Xil {X D ,)) yeD' 

ixAv) y i D ', 

y(£D',n>l 



(39) N^ n (fa = { lx , u {y) y£D',n = l 



For card(I) > 1, 
(40) N^ n {fa) 



n y (e-&>+i»HZ tn (X Dk )) yeD! 
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where Ij = I — {Dj} . Note that we have 

1 



•y 



n Cn (<£i , • • • , <Pk + u, Y 1 , ■ ■ ■ , l Un ) (y)Kx,n {v, dui,..., dv n ) 



where C n = {1, . . . , n}, and the n^. are the operators defined by equations (8) and (9). 



(41) 



Lemma 9. Assume (27). Let 0" = (0i, . . . , fa-i, 4>k + u). 
For n = 1 and y € Dj, 



For n > 2 and y € Dj, 

n—l „ 

N^ n m) = E / n 



m=l 



(42) 



u 

I.v.n 



TD, 



+ii y N^ I )M^ i 'm. 



2Nl^ m ( ( p I )N^'^- m ( ( pjW t I ^mdt ) K x (v,du') 



Proof. If n = 1 then H^i(n) = e~^' u ^ {fi, j x ,u) an d the result is an immediate consequence of 
the basic Palm formula (5)(for card(I) = 1) and the extended Palm formula (11) (for card(I) > 1). 
If n > 2, by the recursive moment formulae (6), 

- t [ N y (e-^*>e-^^n j (X Dfe ,7^))K a;in (i/;^ 1 ,...,^ n _i) 



A + B 



where 



.4 



2-n! 



E n , 



McN 
0,N=£M 



TD, 



/i, . . . ,du n ) 



•%( e -^^)n i6M (^,7^)) 
•%(e^ <X/ '^ > n^ M (^D',7x,, l ))^)^x,n(^;^i 

There are ( n ) possible choices of M in the above expression, with cardinality m. Therefore, by 
rearranging the indices, and using Lemma 6 with r = 2 we get 

n-l 



m=l 



A = E / K x(^du')TI y 



TD, 



%(e _<X/ '^ > n^r(X Dfc , 7 ^)) K x ,„_ m (z, - u'; du x , . . . , dv n - m -i)] dt) ■ 



. (n — m)\ 

The term £> is for card(J) = 1 and for card(I) > 1 is found using the extended moment formula 
(12) as 
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□ 



Now set T x>u = 2 J 'JxmHx^-u' K x (w,dv'). 

Theorem 10. Assume (27). The function j X:l/ is C 4u -superharmonic, and hence lower-semi- 
continuous in y. For R x -almost all v € M c it is in fact an C iu -potential, and satisfies 



(43) 



G 



K x (v;du') 



Proof. Let be a sequence of domains exhausting D. Since Vo k u = u, we have that J\ff k (u) 
e -/oM6)<fc for t < Thug 



71=1 

OO 71— 1 n 



n=2 m=l 



T Di. 



2N t D * («)N et (<J% (H^J dt) K x (u; dv>) 



ww + EE / n 

m=l j=l 1 



4« / 
2/ 1 



n y (H^(X Dk )) + 2 / / 7xjI/ , (dt)l x ,»-v> (6) dt) ^(y; di/) 



n y {H^{X Dk ))+K 



Au i 



The first term is £ 4u -harmonic on by Lemma 9, and the second term is an £ 4u -potential, so 
7xii , is £ -superharmonic on each D^. Thus it is so on D as well. 
Moreover 



/ 



n y (H^(X Dk ))R x (du) = n y (e-( x ^ J (X Dk , 7x , u )R x (dv)) 

= N y ( e -< X *V«> J J X Dk {dw) 1XtU (w)R x {du)) 
= N y (e-< XD *' U > [ X Dk (dw)N w (X D ^0)) 



= e"(f)p y (e-< XD *'«>(X Dfc ,ti)), 

where the last equality follows from the recursive moment formula (7). As we argued in the proof 
of Lemma 5, the stochastic boundary value of u is or oo Py-a.s., therefore P y {e~ {X °^ u ) {X Dk ,«» 
converge to as k — > oo by the dominated convergence theorem. Thus 



as k — > oo. By Fatou's lemma, 
(44) 



N y (H^(X Dk ))R x {dv) ->■ 



liminfN,(^(X Ofe )) = 

k— yoo 
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for R x -a.e. v. But the second term in our expression for j xv (y) 1S monotone in k, therefore 
linifc-^tx) Ny(-ff^ (XD k )) must exist and thus equal to for R x -a.e. v by (44). So, we get that 
7x,u{y) = n^(/ rDfc r X)1/ (£ t ) dt) for R x -a.e. u, for every y. 

Choosing a countable dense set y\, y2, ■ ■ ■ , we therefore have a set Mo D M such that R X (M§) = 
0, and the above equality holds for every j/j and for every v € (A/"o) c . Since both functions are 
lower-semi-continuous, and agree on a countable dense set, it follows that (43) holds for every 

v e (M ) c . □ 

Suppose now that v € M§ where Mo is the -R x -null set described in Theorem 10. Again we may 
assume that if u G M$ then K x>n (u, dv\, . . . , dv n )-ahxvost all {y\, . . . , u n ) are in (M c ) n . Let N y be 
the excursion measure of a super-process whose spatial motion is killed at rate u. In other words, 
for any D' g D, we have 

N y (F(X D ,)) = N y (e-< x ^F(Xjyfj . 

Let y € D be such that j x ,u(,y) < oo. We define a probability on an auxiliary probability 
space Q where there is a branching diffusion on D, and conditional on this branching diffusion a 
Poisson random measure is generated on the infinite product space MP 13 - . We endow MP 13 - with 
the c7-algebra T c d~ generated by the coordinate maps xjj>, D' <s D (That is xd>(uj) = ojd' f° r 
uj € A4 0,D -). Our goal is to construct an M. 0,D - -valued process X = (Xn')D'imD on such that 
the law of X with respect to Q v y will be the same as the law of the exit measures {Xd'}d'<^d °f 
a SBM. First we describe how the branching diffusion evolves: We start a 7-r^-transform of a £ 4u 
process off at y. Since j x>l , is a potential, this process dies before reaching dD. Say it dies at w. 
Then almost surely T x ^(w) is finite. Because, 

u 1 A 00 

n^' "'"(ir 1B , v (f c )=oo) = - — rr / n?(ir ie , w « t )=oor a . iV (^)i c > t )(« 

lx,u\y) Jo 



which must be equal to since Uy (lp x v ^ =00 T XiV (^ t )^(>t) is or oo for any t, and the assumption 
Jx,u(y) = ^yi^x^iCt^ot) < oo implies Ii^(l Tx!V ^ =00 T XiV (^t)lc>t) = for Lebesgue-almost all t, 
making the right side of this equation equal to 0. Note also T XtU > since v 6 A/"q, and therefore 
K x (v, du') is a strictly positive measure. Now we can choose u' at random with distribution density 

2 

(45) - — jr-rTzyM (w)K x (u;du ). 

J- x,v\W) 

Almost surely, both v 1 and v — v' are in J\ffi, and ^y x ,u(w) and J x y(w) are finite so we may start 
two new processes at w, following j x v and ^y X:U -^ transforms of C 4u . Note that we may repeat 
this process infinitely often. This defines a branching particle system. Let Tj denote the measure- 
valued process putting a unit point mass at the historical paths of each particle alive at time t, (A 
historical path of a given particle at time t is the path describing for any s < t the location of the 
particle or whichever ancestor that is alive at time s). We then create mass uniformly along this 
set of particle paths, which then evolves according to the law N(.). Loosely speaking, we want to 
generate a Poisson random measure with intensity 

poo 

/ 4:T t (dz)N Zt (X G -)dt 
Jo 

and add up the resulting measure- valued processes to form X. But we must be careful to represent 
X£>i using only the portions of T corresponding to particles whose historical paths have never left 
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D'. To formulate this, for each t and historical path z we define the following map X t,z : $7 i— > MP D ~ 

X^luS) , = [ X D>{u) \iT D ,{z)>t 
\ otherwise 

Now we generate a Poisson random measure on MP°- with intensity 

/■oo 

X(A) = 4T t (dz)$ at ((X t >*)- 1 (A))dt 
Jo 

which is now well defined for any A G T c d~, since (X t,z )~ 1 (A) is T z ^_ measurable. (Recall, for 
any c G E, JF^- D _ is defined as the cr-algebra generated by {Xd',D' G 0,c G D'}, which is the 
domain where the measure N c is defined) . Adding up the resulting measure- valued processes gives 
us X = (X D >)r» eC , D -- 

More precisely, the n-dimensional transition operators for X with respect to Qy is given by the 
formula 

Qfto) := Q-( e -^^))=g-( e ~/o °4<T i ,N.(i-e-^ 2 ^))) rf t ) _ 

where I = {D 1 , . . . ,D k } C O, and y G D. 

Note that for y ^ D/, Q I y V {4>i) = Qy ,V (4>ii) f° r some j such that y ^ (since all paths in the 
backbone start from y, Xp Z = 0, almost all z for all t, Q y -almost surely). 

The main result of this section is 

Theorem 11. Assume (27), that v G M§, and that J x , u (y) < oo. Then Ny-law of (XD')yeD'mD 
is the same as the Qy-law of (^£>')y&D'&D- 

Proof. Let y G U x m U 2 <s • • • <s Ui .. . s.t. D = U^C/j. It will suffice to show that Qy V {4>i) = 
Ny V (4>i) for J = {Di, . . . , Z?fc} where Dj C = C/j for a fixed i. 

Let D' = Ui. Define T D (dz) as the random measure on D'-valued paths defined by T D (dz) = 
lim^oo f Y t (dz)lr z T D , gd -| where z TD is the path z stopped at T£>>(z). We will write T D ~ n if the 

support of T D ' consists of exactly n paths (that is, exactly n particles of T exit D'). Let 



Q'y^l) :=Q£( e -^*\T B '~n) 
We will show by induction on n and card(/) that 

(46) N I y ^(4>i) = i,Ay)Ql' v ' n (h)- 

for all y G Dj. The theorem then follows by summing on n. 
First observe that, for card(I) = 1 



(47) i,Av)QT n ^) 

and for card(J) > 1, 



i,v,n/x\ _ / lx,u(y) y<£D',n = l 
y<£D',n>l 



(48) i.AyWy^m = i x Ay)Qv' v ' n v?D j: j^k 
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Comparing these equation to equations (39) and (40) we see that equation (46) holds on dD' . 

Let (£t)t>o be an ^-diffusion under the law U y . Uy U and Uy U, ' Jx '' h ' will denote the laws under 
which £ is respectively an C 4u diffusion, and ^ x ^-transform of an C Au diffusion. In what follows 
£ will represent the first particle of the branching backbone which is by construction following a 
7 Xii/ -transform of an £ -diffusion. We let £ be the lifetime of £. Note first that 

[tAr Dl 

= exp(-4 / K£ s )+f%(e-^'> u > - e'^'^^ds) 
J o 

= exp-4 / [u(£ s )-%(l-e- <x o'^)) +Nes (i_ e -(^^))] ds 
Jo 

If n = 1 then 

= n, (lA& Dl )e" M6) ds e" ^ 4 ^ (i-e-<^>)*ggi ( 07 )) 
(49) = n,( 7 ,(£ TD/ )Q^(^)AA4^(^)) 

The first equation is true, because of the strong Markov property of Tt at the first exit time of 
D[ of its first branch £. The second and third equations follow respectively from the definition 
of 7^-transform, and L 4u diffusion. If card(J) = 1, equations (49) and (41) and the fact that 
Ny' u,1 ((f)i) = r y x ,u(y)Qy 1J ' 1 (4>l) on the boundary of Dj = D' implies Ny ,u,1 (4>j) = r y x ,v(y)Qy U ' 1 ($1) 
holds for y G Dj. If we assume Ny' u,1 ((/)i) = j x ,i/{y)Qy U,l {4'i) holds for y € Dj when card(J) = k— 1, 
then equations (48) and (40) and (49) and (41) imply that N y ' u,1 ((f)j) = r ) x ,u{y)Qy U ' 1 (<Pi) for y € Dj 
when card(J) = k as well. 
If n > 2 then 

lx MQy^ Xl ^\T D ' ~ n ) = A + B 

where 

A = -i x ,u{y)Qy{e~ {Xu4n \T D ' ~ n, the first branch of T dies inside D r ), 
B = J x ,u{y)Qy(e-( Xl '* l} , T D ' ~ n, the first branch of T exits Dj). 
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We compute A as follows: 

(fc) 

m=l ' 



n . 

^4 = T^Cy^'^lf^e-^^^^^E/ 



iV(£c) 

■g^e-^'^.T 23 ' ~m)Q^'(e-^'^,T 13 ' ~ n - m) K x {y\ dz/) 



u y ( I r,,te)i,< ? e ^ / f - m 



m=l 



• Qg(e-<*"*'\ T D ' ~ m)Q^'(e-^ A) , T D ' ~ n - m) Jf*(i/; dz/) dt) 

= II,( re-^MW^-Jo'^d-e-^*))^ / 27^(6)7.^(6) 
V ^ m=l ^ 

• Q^(e-<*'^>, T D ' ~ m)Q^ l/ '(e-^>^,T D ' ~ n - m) dz/) di) 

(50) = £ /n,( P'2<^(^)-7, X (6)Q^(e-^),T^'~m) 

m=l ^ ^° 

• T^-^Q^e-** 1 '^, T D ' ~ n - m) dt) ^(i/; di/). 

The first equation is true, because at the lifetime £ of the first particle T branches into two 
new branching diffusions with joint conditional law Q^~ u x Q^~ u given (£c)t<C' an( ^ y i wri ° se 
conditional distribution given (Cc)*<C nas density equal to (45). The second equation follows from 
a well known fact on /i-transforms when h is a potential, see for example Salisbury and Verzani 
(1999) formula 2.2. The third equation follows from the definition of a killed diffusion. 

We compute B as follows: 

B = 7.,,(y)I^^(lc>n, j e-^'^^" <X ' , *' >) *Qjj;(0/) 

(si) = liyM^Q^ydMi^m) 

by first applying the strong Markov property of T at the first exit of Dj of its first branch, and 
then again using the definition of j XjV transform and killed diffusion. 

We have shown previously that for all /, Ny^' l ((f>j) = r )x,y{y)Q I y V ' 1 (4>i) > V £ Dj. Let us assume for 
all /, N y ' u ' m ((pi) = "y x ,u(y)Qy' u ' m (4'i)i V £ Dl-> for all m < n — 1. If card(J) = 1, comparing equations 
(50) and (51) with (42) and using the fact that N y ' u ' n ((j)j) = r y x ,u(y)Qy' U ' n (4>i) on the boundary of 
Dj = D' , we get that N y ' u ' n ((j)j) = r y x ,v(y)Qy l '' n (<Pi) for all y € Dj. If now in addition we assume 
N y ,1/,n ((j)i) = "fx,u(y)Qy' U ' n (4'l) holds for y £ Dj when card(I) = k — 1, our induction hypothesis 
and equations (48), (40) combined with equations (50), (51) and (42) imply that N y ,u,l {(j)j) = 
lx,u{y)Qy' y ' 1 {4>i) f° r y ^ Dj when card(/) = k as well. Hence equation (46) holds for all / and 
y € Dj, therefore the proof is complete. □ 

Above we described the conditional law N^. Now we move to an arbitrary initial measure [i € M.% 
such that H x (fi) < oo, and so need to handle multiple lines of descent starting from time 0. In other 
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words, we are going to describe the distribution of (-XD')A»e.M D /,D'!=D) Pu- Let P^ be a probability 
denned on an axillary probability space under which a random cluster of points (xj,^)™ =1 in 
D x V* is generated according to a distribution proportional to 

(52) —K xn (u; dvi,..., dv n )^ xvi (x\) • • • 7 Xjl , n (x n )/i(dxi) • • • (i(dx n )- 

Note that this makes sense since 

/ K x ^ n (v]dvi,...,dv n )^ x ^ 1 {xi)...^ x ^ n (x n )^{dx 1 )...ii{dx n ) = H^(fi) < oo 

n=1 n - Jxi,...,x n Ji/i,...,v n 

Once the random cluster is generated, corresponding to each point (yi,Ui) in the cluster, a mea- 
sure valued process X 1 begins to evolve following a Qy* law independent of everything else. This 
is consistent with our construction of the law Qy\ since almost surely, each (yi,Vi) will satisfy 
lx,uXyi) < oo and z/j G A/q. In addition, independent from all this, another measure valued pro- 
cess X° evolves following SBM law with spatial motion killed at rate u and initial measure jj,. Let 

x = zr=i xi + x °- 

Theorem 12. Assume (27). Let H x (fi) < oo. Then P^-law of (Xd')^m d i,d'i&d is the same 
as P^-law of (X^ftzMn^D'mD- 

Proof. Let P^ and P^' v denote the transition operators of X and X. That is P^ ,u ((j)j) = 
P M (e-<**.*'>), where I = (D lf . . . , D k ), X x = (X Dl ,. . . , X Dk ), fa = (fa, • • • , fa). P^(fa) is defined 
similarly. 

Let /x be compactly supported in U\ <g U2 <e • • • <£= Ui ■ ■ ■ s.t. D = U^t/j. It will suffice to show 
that Pj l l ' l/ ((j)j) = Pj l ,v ((j)i) for / = {D\, . . . , D^} where Dj € D^. = U{ for a fixed i. 
Recall = H u (fiy 1 P f ,(e~^ Xl ^H u (X Dk )). Note 

^(e-^'^pM) 

/ 00 /• -(X Dk ,u) \ 

= P„ U-^'^^J (X Dk , lx , Ul )---(X Dk , 1XyUn ) \ K Xin {v;dv ll ...,dv n ) 

= Yl - / P » (e^ XD ^(X Dk , 1XtUl ) ■ ■ ■ (X Dk , lx , Vn )J K x , n (u;dv 1 , . . .,dv n ) 

n=l H ' ^ 

Using the extended moment formula with l I,u (x) = 4N X (1— exp — ((Xu 1 , fa) + • • • + (XD k ,(f>k + u))), 
we rewrite the right side as 

00 1 r 

n=l n ' ^ 7r(n) 
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By Lemma 6 we expand the above expression as 



oo 



e 

n=l '"' 7r(n) ' 

K Xtr (u;dui,...,du r ) 

I K x Ju;du u ...,du r ) — 



-^E^E/^;^ 



ni! • • • n r !r! 

n=l r=l ni,...,n r 



) 



1=1 ^ 

r=l r ' t=lrii=l n *' 

oo „ r 

E~ / K x ^u-d Vu ...,dv r )\[^, lx ^Q^{e-^^)) 



oo 

il,u\ 



r=l " i=l 

Note that e - ^'' 7 '"^ is the transition operator of a SBM whose spatial motion killed at rate u, and 
the rest of the expression is the transition operator of T,f =1 Xi, where each Xi evolves according to 
Qy*, where the random cluster of points (yi, z^)?=i is selected according to the density (52). Hence 
the proof is complete. □ 

Note that the extended X-harmonic functions H x were defined by Dynkin in Dynkin (2006b), 
and we have followed this approach throughout. The results of this section should be viewed as an 
attempt to clarify the structure of these extended X-harmonic functions, as well as the structure 
of the conditioned superprocesses that are obtained from them. 
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